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Abstract 

The thermo-economic optimization of an irreversible Stirling heat pump cycle with a detail parametric study for the finite heat capacity 
of external reservoirs is presented in this article. The external irreversibility is due to finite temperature difference between the working fluid 
and the external reservoirs while the internal irreversibility is due to the regenerative heat loss. The Thermo-economic function is defined 
as the heating load divided by the total cost of the system along with the running cost. The Thermo-economic function is optimized with 
respect to the working fluid temperatures and the values for various parameters at the optimal operating condition are calculated. The effects 
of different operating parameters on the performance of the cycle have been studied. It is found that the effect of regenerative effectiveness 
and the economic parameter are more pronounced than that of the other parameters. 

© 2003 Elsevier SAS. All rights reserved. 
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1. Introduction 

The Stirling cycle is one of the important models of re¬ 
frigeration systems for the production of desirable tempera¬ 
tures and the new developments in the design were proposed 
after the new concept of finite time thermodynamics [1,2] 
came into existence. It was remarked that the efficiency of an 
engine operating at maximum power is given by a different 
formula [i] m — 1 — ^/TlJTr], which is always smaller than 
the well known Carnot formula [rjc = 1 — T\J 7’h], where 
7 l and T\\ are, respectively, the temperatures of sink and 
source for heat engines while it is reverse for heat pump cy¬ 
cles] and agrees much better with the measured efficiencies 
of operating installations. Leff and Teeters [3] have noted 
that the straight forwards C-A calculations will not work for 
a reversed Carnot cycle because there is no “Natural Max¬ 
imum”. Blanchard [4] has applied the Lagrangian method 
of undetermined multiplier to find out the COP of an en- 
doreversible Carnot heat pump operated at minimum power 
input for a given heating load. 
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In recent years the performance of the different heat 
engines and refrigeration systems have been investigated 
by a number of researchers using the concept of finite 
time thermodynamics [5-15], while others have applied the 
ecological [16-23] and thermo-economic [24-29] approach 
for various operating conditions. 

In this paper we will investigate the thermo-economic 
optimization along with a detail parametric study of an 
irreversible Stirling heat pump cycle for different operating 
conditions using the concept of finite time thermodynamics. 


2. System description 

It is well known that the working substance of the 
Stirling cycle may be a gas, a magnetic material etc. and for 
different working fluids, this cycle has different performance 
characteristics. When the working substance of the cycle 
is the prefect/ideal gas the Stirling cycle consists of two 
isothermal and two isochoric processes. The schematic and 
T-S diagrams are shown in Fig. 1. This cycle approximates 
the expansion stroke of a real cycle as an isothermal 
process 1-2 with an irreversible isothermal heat addition at 
temperature T c from a heat source of finite heat capacity 
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Nomenclature 


A 

C 

Cf 

COP 

p 

Q 

R 

Ro 

T 

t 

U 


area. 

heat capacitance rate. 

molar specific heat. 

coefficient of performance 

power. 

heat. 


. kW-KC 1 

kJkmol _1 -K _l 


kW 

kJ 


heating load 
gas constant 
temperature 
time. 


. kW 

kJkg _1 -K _1 

. K 

. s 


overall heat transfer coefficient.. kW-m -K 1 


1,2, 3,4 state points 


Greek letters 
e effectiveness 

A volume ratio 

Subscripts 

c cold side/source side 

f fluid 

h, H sink side/heating 

L heat source 

max /m maximum/optimum condition 

R regenerator 



s 


Fig. 1. Schematic and T-S diagrams of an Irreversible Stirling Heat Pump. 

whose temperature varies from 7 li to 7 l 2- The heat addition 
to the working fluid from the regenerator is modeled as an 
isochoric process 2-3. The compression stroke is modeled as 
an isothermal process 3-4 with an irreversible heat rejection 
at temperature Tj, to the heat sink of finite heat capacity 
whose temperature varies from T\\ \ to T\\i. Finally, the 
heat rejection from the working fluid to the regenerator is 
modeled as an isochoric process 4-1, thereby completing the 
cycle. 


As mentioned earlier, the heat transfer processes 1-2 and 
3-4 in a real cycle must occur in finite time. This requires 
that these heat processes must be proceeded through a 
finite temperature difference and therefore, defined as being 
externally irreversible. There is also some heat loss through 
the regenerator, as an ideal regeneration requires an infinite 
regeneration time or area, which is not the case in practice. 

3. Thermodynamic analysis 

Let <2c is the amount of heat absorbed from the source 
at temperature T c and <2h is the amount of heat released to 
the sink at temperature 7h during two isothermal processes, 
then: 

Qh — ThAS — Ch(7h? - ThiRh = CnfiHlTh - ThiRh (1) 
Qc = T C AS — Cl(7li - 7L2)fL = Cl£l(Tli - T c )!l (2) 
where 

AS — nRolnX (3) 

A is the volume ratio, n is the number of mole for the 
working fluid, and Rq is the universal gas constant. Ch, Cl 
and fi. are the heat capacitance rates of sink/source 
reservoirs and heat rejection/addition times and £h and £l 
are the effectiveness of the heat exchangers for hot- and cold- 
side, respectively. Also from the heat transfer theory the heat 
transfers Q h and Q c will be proportional to the Log Mean 
Temperature Difference (LMTD), i.e., 

Qh = C H AH(LMTD) H rH (4) 

Qc = C L A L (LMTD) L r L (5) 

where UhAh and LjAi, are the overall heat transfer 
coefficient-area products and (LMTD)h and (LMTD)l are 
the Log Mean Temperature Difference on sink and source 
side, respectively, and defined as: 

(LMTD)h = [ l7; ' /t ' 1 ' I 7 ' 7 "-' 1 (6a) 
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(LMTD) l = 


(Tli ~ T c ) - (T L2 - T c ) 

In(fel) 


From Eqs. (l)-(6), the effectiveness of the hot- and cold- 
side heat exchangers are obtained as: 


e H = 1 


£L = 1 


where N\ [ and /Vi, are the number of heat transfer units based 
on minimum thermal capacitance rates and are given as: 

( aj C/ h Ah 

V H = ~ciT' Nl = ~cT) 

As these cycles, in general, do not possess the condition 
of perfect regeneration, so it is assumed reasonably that the 
regenerative loss per cycle is proportional to the temperature 
difference of the two isothermal processes as assumed by 
earlier workers [7-10,12-15], 

A<2r = MC f (l — £r) (7], — T c ) (7b) 

where cf is the molar specific heat capacity of the working 
fluid and £r is the effectiveness of the regenerator, which is 
defined as: 


Q regen, actual 
(2regen, ideal 

Ty~T c 
7|, - T c 1 


_ (?4X 
241 
Nr_ 
+ Nr 


T h -T x 
7h - 7’c 


where Nr is the number of transfer unit of the regenerator 
and defined as: 

„ (UA) r 


here Cf is the heat capacitance rate of the working fluid. 

Owing to the influence of irreversibility of finite heat 
transfer, the regenerative time should be finite as compared 
to the two isothermal processes and is given by [5,7-10,12- 
15]: 


fit = U + ?4 = 2a(7h — 7’c) 


p (2h-2l) (2h-2c) 

fcycle (fit + th + ?r) 

(13) 

„ 2h 2h 

AR = - = - 

Uycle (fit + *L + r R) 

(14) 

7?h 2h 

COP h = = 

p (2h-2c) 

(15) 

Using Eqs. (11)—(16) we obtain: 


2c 2h q 

(16) 

Thus, from Eqs. (11)—(16), we have: 


(x — 1) 

P = --—p-- 

U H (*y-r H1 ) + *L(r L1 -y) + b ^ x 1)] 

(17) 

[x — ai(x — 1)] 

AH = -;-j- 

UhCvv-Thi) + k L (T Ll -y) + b ^ X 

(18) 

[x — ai(x — 1)] 

COP H = ., 

(x - 1) 

(19) 


where ku = £hC h , k L = £ L C L , x = T h /T c , y = T c , b\ = 
2a/(«7?olnA) and a\ — cf(l — £r)/7?o lnA. 

The objective function of thermo-economic optimization 
proposed by earlier workers [24,25] is given by: 


F = 


2h 

Ci + Ce 


( 20 ) 


where C; and C e refer to annual investment and energy 
consumption costs. 

The investment cost was considered the costs of the 
main system components that are the heat exchangers 
and the compression and expansion devices together. The 
investment cost of the heat exchangers is assumed to be 
proportional to the total heat transfer area [24,25], On the 
other hand, the investment cost due to the compression and 
expansion devices is assumed to be proportional to their 
compression/expansion capacities or the required power 
input. Thus the annual investment cost of the system can be 
given by: 


where a is the proportionality constant, which is indepen¬ 
dent of the temperature difference but depends on the prop¬ 
erty of the regenerative material. Thus, the total cycle time 
Cycle IS: 

Cycle = OH + th + fit) (10) 

When the irreversibilities mentioned above are taken into 
account, the net amounts of heat released to the sink and 
absorbed from the source are: 

2h = 2ii-A2r (11) 

2l = 2c ~ A2 r (12) 

It is well known that the power input, heating load 
and coefficient of performance (COP) are the important 
parameters of the heat pump. Using the above equations we 
have: 


C; = a(An + Al + Ar) + b2P 

— o(Ah + Al + Ar) + (>2(2h — 2c)/ Cycle (21) 

where the proportionality constant for the investment cost 
of the heat exchanger, a, is equal to the capital recovery 
factor times investment cost per unit heat exchanger area 
and the proportionality constant for the investment cost for 
the compression and expansion devices, bi, is equal to the 
capital recovery factor times investment cost per unit power 
input. The initial investment cost is converted to equivalent 
yearly payment using capital recovery factor [24,25], The 
annual energy consumption cost is proportional to the power 
input, i.e., 

C e =b^P = 7>3(2h - 2c)Acycle (22) 

where the coefficient, bj, is equal to the equivalent annual 
operation hours corresponding to minimum power input 
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times price per unit energy [24,25]. Substituting Eqs. (21), 
(22) into Eq. (20) we have: 


F = 


<2h 


a (Ah + Al + Ar) + b(Q h — Qc)Acycie 
where b — bi + b^,. Thus from Eqs. (21)-(23) we have: 

[x — a i (x — 1)] 


(23) 


bF = 


k\ki[ 


kn(xy-T m ) 


k L (T u -y ) + b t (*-!)]+ (*-!) 

(24) 


where k\ — a/b and 
CffiR 


h = 


Ur(1 - £r) 


ln(l - £h) - jj- ln(l - £l) 
Uh U l 


It can be seen from Eq. (24) that bF is a function of 
two variables x and y [as other parameters are constant for 
a typical set of operating conditions]. Thus optimizing bF 
with respect to ‘ y’ viz. — 0 yields: 


x( T{ \ -y) = . 


(xy - Phi) 


Substituting Eq. (25) into Eq. (24) we have: 


bF = 


[x — a i(x — 1)] 
[(.rTu-Vni) +b ^ x ~ 


(25) 


(26) 


where 



ks — k i A '2 A: 4 


k4 = a+ k 3 ) 

\kn k\k^) 
and Z ?4 = 1 + kik 2 b\. 


Optimizing 'bF' with respect to x viz. klL = 0 yields: 


*opt= {^PhiPli + [7hi7li[6 4 7hiA- 5(1 -fli) 

+ b 4 aiT Ll k 5 -aikj(l — «i)]* /_ } 

x {T Ll [b 4 T Ll - k 5 (l - (27) 

The maximum objective function and the corresponding 
power input, heating load and COP can be expressed as: 


UI7 l*opt ^l(-^Opt 1)J 

b^max — u 

+**<*». PI 

(28) 

p C^opt — 1) 

[<,„,£-*.,>+*■<** pi 

(29) 

D l*opt a \ C^Opt 1)J 

^H,m — k 

1)] 

(30) 

„„„ U'opt — fl l(^opt — 1)J 

COP H ,m — 

l-^opt 1) 

(31) 


3.1. Special cases 


(1) When £r = 1.0, i.e., the cycle possesses the condition 
of perfect regeneration, the maximum thermo-economic 
function and the corresponding power input, heating load 
and COP are given by: 


bFmw = ■ 


•^opt 


Aopt^5 


Pm — 


Pllin = 


[(^&I)+^ opt-D] 
C^Opt 1) 

f (-Voptru-Tm) + b ^ x opt - 1)] 

•^Opt 


^ t v op . /Ti* 5 y„i) + /; i(- v -.pt i)] 


COPn.m = 


•^opt 


C^Opt 1) 


where 


•*opt : 


7hi b 4 T L1 + V^sPli 
Pli [^ 4 Pli - ^ 5 ] 


(32) 

(33) 

(34) 

(35) 

(36) 


In this case not only the maximum thermo-economic func¬ 
tion but also the corresponding heating load of the Stirling 
cycle is similar to that of an equivalent Carnot cycle. How¬ 
ever, for finite regeneration time/area £r should be less than 
unity. It shows that in the investigation of the Stirling cy¬ 
cle it would be difficult in obtaining new conclusions if the 
regenerator losses were not considered in the analysis. 

(2) If /r = 0, i.e., the time of the regenerative processes 
is negligible as compared with that of the two isothermal 
processes, the maximum thermo-economic function and the 
corresponding power input and heating load are: 


bPmax = 


-*opt <41 (VdpL 1) 
Q-vopiTH-rHi) + ^ 


Pm — 


(X, 


Opt ' 


i)(* 0 p t r li 


Phi) 


P|[m = 
where 


CvTipl 

L*opt tt 1 (vdpL 1) J (xopt Pjj Phi) 
CvPopt 


(37) 

(38) 

(39) 


Jfopt = {PhiPli + VPhiPli [PhiMI - at) 

+ oi7li^5 -aikj(l -ai)] 1/2 } 

x{P L1 [P L1 -^5(l-a 1 )]r l (40) 

From the above results, we conclude that the Stirling cycle 
has different performance characteristics for different con¬ 
figurations. Thus, the expressions given by Eqs. (17)-(20) 
are more general and include all the major irreversibilities 
associated with this cycle. 

(3) When Ch —*■ 00 and Cl -* 00 , i.e., the heat capaci¬ 
tance rates of external the reservoirs are infinite, the Stirling 
heat pump is operating between two infinite heat reservoirs. 
In such a case, the maximum thermo-economic function and 
the corresponding power input and heating load are: 
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bFmm — 


Pm — 


■'-"opt 0 1 (■'-Opt 1) 

( (3^7^55) + ^ 4 ^ opt “ ^ 


(x, 


Opt ‘ 


1) 


( (JCppt7L -TH) +b4 ^°P l 1 ^ ^ 


%n = 


L*opt a\ Oopt 1 )J 


&5-*opt 


(^opt7L 7 h) 


-b 4 (x opt - 1)) 


where 



Jfopt = {b^TuT-L + [7h7l[&47h^ 5(1 - «i) 

+ ^ 4 ai7 l& 5 - aifcjd - fli)]]'^} 
x {T L [b 4 T L - k 5 (l - ai)]}~ 1 (44) 

while 

% = ,/^S and t 4 =('-d = + -J = y 
V VV^hAh JUlA l/ 

(4) When 7 r = y (?h + <l), he., the time of regeneration 
is directly proportional to the time of two isothermal 
processes, the maximum thermo-economic function and the 
corresponding power input and heating load are: 

(45) 

(46) 

(47) 

where 


bF mm = 


-V'opt O I (A'opt 1) 

(l+K)A- 5 -V op , _ 

M-VopiTLl-Thl) + '■oPt ^ 


Pm — 


P\\m = 


(■'Opt 1 )(-'opt'/’u f ill) 

(1 + y)A: 5 x 0 pt 

L-'opt a I (-'opt 1)J (■'opt7’[,i 7) 11 ) 


(1 + y)fcs*opt 


■t-'opt = {ThiTli + [7hi7li[7hi^ 6(1 + y)(l - at) 

+ «i7li* 6 (1 + y) - + y)(1 - «i)]]‘ /2 } 

x {i[! —& 5 (1 + y)(1 — fl i)]} _1 (48) 


4. Discussion of results 


(a) 



(b) 


Fig. 2. (a) Heating load and power unput vs Ic\. (b) Objective function and 
COP vs k\. 


4.1. Effect of k\ 

The effects of economical parameter (k i = a/b) on 
the heating load, power input, COP and thermo-economic 
objective function are shown in Figs. 2(a), (b). It is seen 
from these figures that the COP increases while the power 
input and the objective function decrease whereas the change 
in the heating load is very small as k\ increases. Since our 
prime object is to have higher thermo-economic function for 
a typical set of operating condition. Hence, it is better to 
have lower value of k\ for the maximum objective function 
from the point of view of thermodynamics as well as from 
the point of view of economics. Again, the effect of k i is 
more pronounced for the maximum objective function and 
less pronounced for the heating load. 


In order to have the numerical appreciation of the results 
for the Stirling heat pump, we continue to investigate 
the effects of the heat sink and heat source temperatures 
(7j[i, 7 li), the effectiveness of the heat exchangers (fin, ep 
and fi k ), the economic parameter (k\) and heat capacitance 
rates (Ch and Cl). While the effect of each one of these 
parameters is examined and the rest of the parameters are 
kept constant as (k i = 0.50, b\ = 0.01, Ch = Cl = Cr = 
1.0, £h — £l — £r — 0.75, Thi = 330 K, 7jj = 290 K, 
X — 2.0, Ch = Cl = 1.0 kW/K) and the results obtained are 
as follows: 


4.2. Effect of effectiveness 

The variation of ( bF ) with the COP for different effec¬ 
tiveness is shown in Fig. 3(a). It is seen that the objective 
function first increases and then decreases as the COP in¬ 
creases. Thus there is an optimum value of COP at which 
the objective function attains its maximum. So it would be 
better to choose the point of maximum objective function in 
the design of the Stirling heat pump from the point of view 
of thermodynamics as well as from the point of view of eco¬ 
nomics. 
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COP m 


(a) 



Effectiveness 


(b) 



Effectiveness 

(C) 

Fig. 3. (a) Objective function vs COP. (b) Heating load vs effectiveness, 
(c) Power unput vs effectiveness. 

Figs. 3(b), (c) show the effects of hot- and cold-side and 
regenerative effectiveness on the heating load and power 
input for a given set of operating conditions. It is seen from 
these figures that both the parameters, i.e., the heating load 
and the power input increase, as the effectiveness on the 
either side heat exchangers increase. But the change in the 
heating load and the power input with the hot- and cold-side 
effectiveness is very small as can be seen from these graphs. 
So it would be better to select the optimum condition from 
the point of thermodynamics as well as from the point of 
view of economics. Again, the effect of regenerative side 
effectiveness is more pronounced than the hot- and cold-side 
effectiveness on all the parameters. Also the effects of hot- 
and cold-side effectiveness are almost the same for these 
parameters, so the two curves overlap as can be seen from 
these figures. 



Thi (K> 


(a) 



(b) 


Fig. 4. (a) Heating load and power unput vs 7j]i ■ (b) Objective function and 
COP vs r H1 . 

4.3. Effect of sink and source temperatures 

Figs. 4, 5 show the effect of sink and source temperatures 
on the heating load, COP, power input and the objective 
function. It is seen from these figures that as the sink side 
temperature increases, the heating load increases while the 
COP and the objective function decrease. On the other hand, 
the heating load decreases while the COP and the objective 
function increase, as the inlet temperature on the source 
side reservoir increases. Also there is a small change in the 
power input in both the cases. Since our prime object is to 
get higher performance with the lower investment, so it is 
better to have higher value of 7 li and lower value of Thi . 
Again, the effects of Thi and 7 li are more pronounced for 
the objective function and less pronounced for the heating 
load. 

From the above results we conclude that the sink side 
inlet temperature should be lower and the source side inlet 
temperature should be higher for better performance of 
these cycles. But it is little bit difficult to decrease the 
sink side temperature rather than to increase the source side 
inlet temperature and can be done by using low-grade-heat- 
source such as solar energy [30] or the waste heat [31] from 
industrial/power plant. 

4.4. Effect of heat capacitance rates 

The effects of heat capacitance rates (Ch and Cl) on the 
heating load, power input, COP and the thermo-economic 
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Fig. 5. (a) Heating load and power unput vs 7 li . (b) Objective function and 
COP vs 7 li . 




C H or C L (kW/K) 

(b) 

Fig. 6. (a) Heating load and power unput vs Ch or Cl- (b) Objective 
function and COP vs Ch or Cl- 


function ( bF ) are shown in Figs. 6(a), (b). It is seen from 
these figures that the COP increases while the objective 
function decreases as the heat capacitance rate on either side 
heat exchanger increases. On the other hand there is a small 
change in the power input and heating load, yet the graphs 
show the trends. Again the effect of Ch and Cl are equal on 
all the parameters. 

Since, the higher value of heat capacitance rate forces 
the working fluid to absorb/reject heat at lower temperature 
by decreasing the external irreversibility. Thus, it would 
be better to have the higher values of heat capacitance 
rates on both side heat exchangers for better performance 
from point of view of thermodynamics while it is better 
to have the lower values of Ch and Cl from economic 
point of view, which is the prime object of this analysis. 
Hence, the optimum operating condition should be chosen 
carefully for better performance from the point of view 
of thermodynamics as well as from the point of view of 
economics. 


5. Conclusions 

The Thermo-economic optimization along with the de¬ 
tailed parametric study of an irreversible Stirling heat pump 
cycle has been evaluated for the finite heat capacity of ex¬ 
ternal reservoirs. It is found that there is an optimal value of 
COP for a typical set of operating conditions at which the ob¬ 
jective function attains its maximum. The thermo-economic 
function is found to be increasing function of the source side 
inlet temperature and the regenerative effectiveness while 
it is found to be decreasing function of the economic pa¬ 
rameter, sink side inlet temperature and heat capacitance 
rates. 

The effects of regenerative side effectiveness and the 
economic parameter are found to be more than those of the 
other parameters from the point of view of thermodynamics 
as well as from the point of view of economics. Moreover, 
the economic parameter, the source side inlet temperature 
and the heat capacitance rates can be used as the important 
criteria in the design of an irreversible Stirling heat pump 
cycle and other similar cycles too. Again, it is also found 
that the inlet temperature on the source side reservoir should 
be little bit higher for better performance of the cycle. It can 
be done by using low grade energy like solar energy and/or 
the waste heat from the industrial and power plants for the 
overall better performance of this cycle. 
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